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Abstract. An interesting accident of nature is that the peak of the cosmic ray 
spectrum, for both protons and heavier nuclei, occurs near the pion production 
threshold. The Boltzmann transport equation contains a term which is the cosmic 
ray flux multiplied by the cross section. Therefore when considering pion and kaon 
production from proton-proton reactions, small cross sections at low energy can be as 
important as larger cross sections at higher energy. This is also true for subthreshold 
kaon production in nuclear collisions, but not for subthreshold pion production. 
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1. Introduction 



The existence of the pion was predicted by Yukawa [lj in 1935 and discovered in 
the comsic ray spectrum by Powell [2] in 1947. Seventy years later one would have 
thought that of all reactions, pion production in proton proton collisions was thoroughly 
understood. There has been considerable interest recently in meson production from 
nucleon nucleon collisions near the meson production threshold j3l HI [5], [6]. Pion 
production has been studied in the most detail. The reason these studies raised 
such interest was because the first calculations of total cross sections disagreed with 
experiment by a factor of five [5] and the differential cross sections for spin observables 
had the wrong shape, being convex instead of concave [7j. 
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Particle reaction transport codes, such as GEANT [8] and FLUKA [5] are widely 
used in the design of accelerator experiments and the simulation of particle detectors. 
While typical cross sections away from threshold reach tens of millibarn, the cross 
sections near threshold are about a thousand times smaller, typically in the microbarn 
region. The fact that theory and experiment disagreed so much in the threshold region 
was therefore of no real concern for particle simulation codes because the cross sections 
in the threshold region are so small. 

Transport codes are also widely used in studies of cosmic rays. If one knows the 
cosmic ray spectrum incident on top of the Earth's atmosphere, then one can deduce the 
cosmic ray spectrum observed on the ground by transporting through the atmosphere. 
Similarly if one knows the cosmic ray spectrum incident on a spacecraft wall then one 
can deduce the radiation environment inside a spacecraft [TUj. In cosmic ray research 
one turns this around to deduce the incident spectrum from knowledge of the ground 
based spectrum. 

Transport codes are concerned with solving the Boltzmann transport equation. In 
the continuous slowing down and straight-ahead approximations, the one dimensional 
Boltzmann equation can be expressed [10] as 

4>j{x,E) = f dE' ^2 T,j k (E, E')cj) k (x, E') (1) 

J k 

where S(E) is the stopping power, £ is the macroscopic cross section and <fi is the flux of 
particles at a given depth x and energy E. The macroscopic cross section is related to 
the microscopic cross section a by X = pa, where p is the target number density. Note 
that the flux is multiplied by the cross section a. In a particle accelerator, the incident 
flux is normally a beam (constant flux) of particles with a well defined energy. In cosmic 
ray physics the incident flux is represented by the incident cosmic ray spectrum, which 
is discussed in the review by Simpson [TT] . In Figure 5a of that review one can find the 
incident spectrum for protons and other nuclei. The peak of the spectrum occurs at an 
energy of about 300 MeV. The peak is caused by solar modulation of the interplanetary 
magnetic field and is known to move around somewhat depending on solar activity. 
Note that the peak of the spectrum occurs quite near the pion production threshold of 
290 MeV. This is an interesting accident of nature. Above this energy the spectrum 
starts falling. 

Now the Boltzmann equation contains flux multiplied by cross section. Therefore 
just because a cross section is small does not mean that it is unimportant. A small cross 
section might get multiplied by a large flux and a large cross section might get multiplied 
by a small flux, leading to approximately equal contributions from a transport point of 
view. This is possible when the incident flux is the cosmic ray spectrum at the top of 
the Earth's atmosphere. 



-- m S A E) + E A E) 
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2. Pion and kaon production in proton-proton collisions 

Consider the exlusive reaction for production of neutral pions, 

p + p^p + p + n° (2) 

Sample values of cosmic ray proton flux [11] and microscopic cross section are given in 
Table 1, where the values are also multiplied together (we call this importance), as in 
the Boltzmann equation. It can be seen that the small cross sections near threshold 
can be of similar importance to larger cross sections away from threshold, the reason 
being due to the fact that the cosmic ray spectrum has a peak near threshold and falls 
steadily at higher energies. For example, the cross section at 375 MeV is 40 fjb, while 
the cross section at 7 GeV is 1.7 mb, yet they give approximately equal contribution to 
the importance factor. Of course these do not compare to the importance of the cross 
sections at say 700 MeV and 1 GeV. However in a good transport code one will always 
include the cross section at a variety of energies and one would certainly go out to 7 
GeV. The point is that if one does include such higher energy cross sections, then one 
may also need to include the cross sections near threshold as well because they are of 
comparable importance. 

Table 1. Neutral pion production from proton-proton collisions. (Importance = Flux 
x a.) The numbers in square brackets after the cross sections indicate the Reference 
number from where the cross section was taken. The peak in Flux occurs near the pion 
threshold. 



T Flux [TT] a Ref. Importance 

( MeV ) (m 2 s/sMev) (l^) ( m 2 sr^Mev) 



290 «- threshold 




325 


2.0 


350 


2.0 


375 


1.8 


400 


1.8 


700 


1.3 


1000 
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3000 


0.24 


7000 


0.06 


11500 


0.01 



7.7 
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17 
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40 
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Similar statements can be made concerning kaon production as shown in Table 2, 
where the lowest threshold reaction for producing kaons, 

p + p^K + + A + p (3) 

is considered. It can be seen that the importance of near-threshold reactions is compara- 
ble to reactions at all energies. Therefore it may be important to consider near-threshold 
production of kaons. 

Table 2. Kaon production from pp — > K + Ap. 



T Flux jTT] a Ref. Importance 

(^ eV ) (n^s/sMev) ( m 2 sr^Mev) 



1.58 <— threshold 




1.60 


0.6 


1.82 


0.5 


1.90 


0.4 


2.06 


0.3 


2.85 


0.2 


4.83 


0.1 


7.29 


0.06 


30 


.002 


50 


0.0005 



0.16 


[13J 


0.1 


7.4 


m 


3.7 


8.6 


m 


3.4 


16.5 


M 


5.0 


50 


m 


10 


50 


m 


5 


50 


m 


3 


16.6 


m 


0.03 


16.7 


1161 


0.01 



3. Sub-threshold pion and kaon production in nucleus-nucleus collisions 

Based on the above considerations, one might ask about the subthreshold pion pro- 
duction observed in nucleus- nucleus collisions some decades ago [T7]. Here the word 
subthreshold refers to pions produced in nuclear collisions where the nucleus kinetic 
energy is below the nucleon-nucleon threshold. The typical way to express energies in 
nucleus-nucleus collisions is using A MeV. An energy of 290 A MeV would represent 
the nucleon-nucleon threshold. Simpson [TT] also provides the cosmic ray spectrum for 
nuclei, such as He, C, Fe. Again the spectrum peak for these nuclei is right near the 
pion threshold. If the peak were below the pion threshold then sub-threshold pion pro- 
duction in nuclear collisions would be important for the same reasons as given above. 
However because the nuclear peaks remain near the pion threshold then subthreshold 
pion production is not important in the sense described previously. This can be seen 
from Table 3 where total cross sections for neutral pion production in Carbon-Carbon 
collisions below the pion threshold are presented. Due to a lack of total cross section 
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data, the experimental cross sections without Reference numbers next to them were 
taken from Table 1 and multiplied by 144, which is the product of the mass numbers 
of the projectile and target nuclei. Therefore they are approximate only. Again these 
cross sections are multiplied by the Carbon cosmic ray flux and it can be seen that sub- 
threshold cross sections are not important compared to cross sections at higher energies. 
Therefore it is probably safe to neglect sub-threhsold pion production in transport codes. 

Table 3. tt° production from C-C collisions. Note that the peak in Flux occurs near 
the pion threshold. 



T Flux [TT] a Ref. Importance 

(AMeV) ( m2 s^AMev) ( mb ) ( m 2 sr^AMeV' 



60 


5 


X 


10~ 3 


1.7 x 10~ 3 pa 


m\ 


8.5 xl0~ 6 


74 
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10" 3 


8.5 x 10" 3 [TBI 


m\ 


5.1 xlO" 5 
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18.9 x 10" 3 [TBI 
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1.3 xl0~ 4 


290 <- threshold 
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0.04 


400 
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12 




0.07 


700 
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10~ 3 


288 




1.2 
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10~ 3 
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1.7 


3000 


5 
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10" 4 


432 
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7000 
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X 


10" 4 


245 




0.02 


11500 


2 


X 


10~ 5 
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0.003 



It is interesting to look at subthreshold kaon production because the kaon threshold 
is well above the pion threshold. Data are presented in Table 4. Again due to a lack 
of total cross section data, the experimental cross sections without Reference numbers 
next to them were taken from Table 2 and multiplied by 3364, which is the product of 
the mass numbers of the projectile and target nuclei. It can be seen that subthreshold 
cross sections are as important as cross sections near 10 GeV. Therefore one can make 
the interesting conclusion that subthreshold kaon production may need to be included 
in transport codes. A better knowledge of the nucleus-nucleus total cross sections is 
necessary in order to make more definitive statements. 
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Figure 1. The percent flux difference of charged muons and pions in the Mars 
atmosphere [24J at a depth of 38.9 from the 1977 GCR spectrum [23] 

4. An example transport calculation 

We have added this section on the advice of a referee, who wished us to do a sample 
transport calculation by artificially multiplying the cross section near threshold by a 
factor of two to see what the effects might be. The transport calculation was done with 
MESTRN [21], an extension of the NASA transport code HZETRN [22] that includes 
the production of charged pions and muons. One improvement that was made to the 
published version of MESTRN is the pion production spectral distributions from proton- 
proton collisions are now calculated by numerically integrating the Lorentz invariant 
differential cross sections of Badwar et al. [23J. This was done to achieve better 
numerical convergence for low density materials and to allow for transport through 
variable density materials. 

The Martian atmospheric model of De Angelis et al. [23] and the primary galactic 
cosmic ray (GCR) model of Badwar and O'Neill [25] were used as inputs for the transport 
calculation. Two runs of the transport code were performed. For one run, the cross 
section for the production of pions was artificially multiplied by a factor of two from 
threshold to a kinetic energy of 2 GeV. 

Figure 1 shows the percent flux difference of charged pions and muons after 
transport through 38.9 of Martian atmosphere corresponding to an approximate 
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altitude of -10 km. A negative altitude indicates a valley or impact basin. The percent 
flux difference is the flux of the MESTRN run with the cross section for pion production 
enhanced by a factor of two minus the nonenhanced flux. This value is then divided by 
the nonenhanced flux to give the percent variation due to the enhancement of the pion 
production cross section near threshold. 

The enhanced pion production cross sections led to an enhanced flux of charged 
pions and muons by approximately 20-30% in the kinetic energy region up to 1 GeV. 
The enhancement fades quickly after 1 GeV. We might expect this effect to be larger, 
but the produced pions decay quickly into muons which subsequently decay also. Muon 
decay in MESTRN is treated simply as a loss term. Subsequently, some of the muons 
that are produced from the increased pion production near threshold are not accounted 
for in these results. 



Table 4. K + production from Ni-Ni collisions. The approximate total experimental 
cross sections are from Reference [20]. To obtain the flux for Ni, the flux of Fe was used 
but multiplied by 0.05, which is the average ratio of the abundance of Ni to Fe over all 
energy intervals [TT] . 



T Flux a Ref. Importance 

(AGeV) ( m 2 sr s AMeV ) ( mb ) ( m 2 sr^AMeV' 



0.8 1.1 xl0~ 5 

1.0 1 xl0~ 5 
1.58 ^ — threshold 

1.60 1 xl0~ 5 

1.8 5 xlO- 6 

1.82 5xl0~ 6 
1.90 4.5xl0" 6 
2.06 3.5xl0~ 6 
2.85 2.5xl0~ 6 

4.83 lxl0~ 6 
7.29 4xl0~ 7 
30 2.5xl0^ 8 
50 5xl0~ 9 



.9 [20] 1.0 xlO" 5 

2.7 [20] 2.7 xlO" 5 

.5 5 xl0~ 6 

57 [20] 2.9xl0- 4 

25 1.3xl0- 4 

29 1.3xl0~ 4 

56 2.0xl0- 4 

168 4.2xl0~ 4 

168 1.7xl0- 4 

168 6.7xl0- 5 

56 1.4xl0" 6 

56 2.8xl0~ 7 
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5. Conclusions 



Of course particles with higher energy will be more penetrating than those of lower 
energy. This aspect of transport is not considered in the present paper. All that is 
considered is the product of the cross section with flux, which we called importance. 
We have shown that the small meson production cross sections near threshold can be 
just as important as larger cross sections at higher energy when transporting cosmic 
ray particles, because the shape of the cosmic ray spectrum enhances low energy cross 
sections, as used in the Boltzmann transport equation. In proton-proton reactions, 
this is true for pion production and especially true for kaon production. This leads us 
to speculate that near threshold production of heavier hadrons may also need to be 
included in cosmic ray transport codes, because their thresholds are all above the pion 
threshold and the cosmic ray spectrum falls steadily as energy increases. 

On the other hand this paper has shown that subthreshold pion production in 
nucleus-nucleus collisions is not important for cosmic ray transport because the energy 
of these reactions occurs below the peak of the nuclear cosmic ray spectra. However 
sub-threshold kaon production probably may need to be included in transport codes 
and again we speculate that sub-threshold production of heavier hadrons may need to 
be included. 

It is important to note that the present paper does not prove that the reactions 
considered must be included in all cosmic ray transport codes. All we have done is 
to suggest that certain reactions with small cross sections at low energy have similar 
importance factors compared to larger cross sections at higher energy, and therefore, 
under certain circumstances, these might need to be included in transport codes. The 
circumstances will depend on the incident spectrum, the particles one is transporting, 
the medium through which transport is taking place and the particles of interest after 
transport through the medium. 
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